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ABSTRACT 

Folding mechanisms in which secondary structures 
are stabilized through the formation of tertiary inter- 
actions are well documented in protein folding but 
challenge the folding hierarchy normally assumed 
for RNA. However, it is increasingly clear that 
RNA could fold by a similar mechanism. PSabc, 
a small independently folding tertiary domain of 
the Tetrahymena thermophila group I ribozyme, is 
known to fold by a secondary structure rearrange- 
ment involving helix P5c. However, the extent of this 
rearrangement and the precise stage of folding that 
triggers it are unknown. We use experiments and 
simulations to show that the P5c helix switches to 
the native secondary structure late in the folding 
pathway and is directly coupled to the formation 
of tertiary interactions in the A-rich bulge. P5c 
mutations show that the switch in P5c is not rate- 
determining and suggest that non-native inter- 
actions in P5c aid folding rather than impede 
it. Our study illustrates that despite significant 
differences in the building blocks of proteins and 
RNA, there may be common ways in which they 
self-assemble. 

INTRODUCTION 

In RNA folding, the stability gap separating the free 
energy of the native state from other non-native structures 
is often smaU (1-3). Thus, misfolded conformations in 
which even the secondary structure contains errors can 
form with substantial probabihty (4). One way to correct 
errors in the secondary structure is by couphng base 
pair reorganization with RNA tertiary folding, either 
by realigning entire helices (5-7) or by shuffling 



non-Watson-Crick base pairs within a helix (8,9). 
Although it has been known for some time that certain 
base pairs are not stable until the RNA tertiary structure 
forms (10), tertiary interactions may also be used to 
correct non-native secondary structures that are present 
in the RNA. Direct couphng between RNA secondary 
and tertiary structure challenges the conventional view 
(3,11) that native secondary structures always form 
before the estabhshnient of tertiary interactions. 

The ease with which RNAs form alternative stable struc- 
tures raises the possibility that tertiary interaction-driven 
base pair switching, and not merely helix stabilization, 
might be more widespread than commonly assumed. To 
explore this possibihty, we investigated the refolding 
kinetics of the P5abc RNA that switches its secondary struc- 
ture when inter-helix tertiary interactions are formed 
(Figure 1) (6,12). The PSabc RNA is an independently 
folding tertiary domain of the Tetrahymena thermophila 
group I ribozyme (13,14) which is one of the fastest 
regions of the ribozyme tertiary structure to form (15). 
P5abc is composed of three helices (paired regions P5a, 
P5b and P5c) that interact via a rich web of tertiary 
contacts in the folded state (16). One set of tertiary inter- 
actions is formed between an A-rich bulge in helix P5a 
and the loop of hehx P5c and another is formed by 
tandem G-A pairs at the junction between helices P5a, 
P5b and P5c. 

A nuclear magnetic resonance (NMR) structure of a 
truncated form of the PSabc domain (tPSabc) showed 
that the PSabc secondary structure in the absence of 
Mg^^ (Figure lA) differs significantly from that of the 
Mg^^-bound state (Figure IB), implying that tertiary 
folding is coupled to a reorganization of the RNA second- 
ary structure (6). First, the base pairs in PSc change 
register, expanding the stable GNRA tetraloop in the 
extended form of PSc to a less stable pentaloop in the 
native state that positions U168 to interact with G188 in 
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Figure 1. Structures of the tP5abc RNA. (A) The secondary structures (top) and ribbon diagrams (bottom) of the extended tPSabc from NMR (6) 
are compared with (B) the folded P5abc structure from X-ray crystallography (22). Mg^^ ions visible in the X-ray structure are represented by 
spheres. Nucleotides are numbered discontinuously because 5 bp of the natural P5b are missing in tPSabc. Non-Watson-Crick base pairs are 
represented by Leontis/Westhof symbols (49). U185 was replaced with the fluorescent analog 2AP (cyan), which unstacks in the folded RNA. 
Base-pairing free energies are from mfold (35). (C) RNase Tl probing of tPSabc secondary structure at 25 and 50°C, with no Mg''^ (extended) or 
2.5 mM Mg^^ (folded). Lane C, no digestion control; lane T, RNase Tl sequencing ladder. 



P5a (6,12,17). Second, reorganization of the three-helix 
junction replaces base pairs A139-U179 and G164-U177 
with tandem G-A pairs (A139-G164 and A140-G163) in 
the native from of P5b. 

Although the native secondary structure of P5abc is 
less stable than its extended secondary structure, this is 
compensated by the favorable free energy of tertiary 
interactions in the presence of Mg^^ ions (17,18). In 
addition to non-specific charge stabiUzation (19,20), 
Mg^"^ ions specifically coordinate a 'corkscrew' in the 
backbone of the A-rich bulge (magenta spheres; Figure 
IB), while four additional ions interact with P5c and with 
tandem G-A pairs in P5b (orange spheres; Figure IB) 
(21,22). Evidence for specific binding of Mg'^^ to P5abc 
RNA in solution comes from NMR, phosphorothioate 
substitution-interference (12,21) and refolding experi- 
ments (23,24). 

Force denaturation experiments confirmed the import- 
ance of tertiary interactions for global stabihty of the 
P5abc domain, as high forces were required to rupture 
the P5abc structure in Mg'^^ (25). More importantly, 
when the force was relaxed, refolding occurred in two 
steps by a kinetic partitioning mechanism (26). Of particu- 
lar relevance is that the single-molecule experiments 
showed that folding of P5abc, upon release of force, 
involves intermediates regardless of the Mg^^ concentra- 
tion. However, such experiments provide no information 



about the structure of the intermediates nor their location 
along the refolding pathways. 

Herein, we provide evidence that helix switching occurs 
in the same step of folding as the formation of tertiary 
interactions in P5abc. Using a fluorescence assay, we 
show that the solution refolding pathway contains at 
least one intermediate. However, experiments and simula- 
tions show that the base pair shift in P5c does not limit the 
rate of folding. Instead, our results show that switching of 
the P5c helix occurs in the same folding step as the for- 
mation of tertiary interactions in the A-rich bulge and that 
both of these events follow the slow step of folding. Our 
results show that the energy landscape of even the rapidly 
folding P5abc subdomain is rugged, exhibiting many of 
the complexities found in the folding kinetics of larger 
ribozymes. 



MATERIALS AND METHODS 

Coarse-grained dynamics simulations 

Coarse-grained simulations used the three interaction site 
(TIS) model (Supplementary Figure SI) in which each nu- 
cleotide is represented by three interactions sites corres- 
ponding to the centers of masses of the base, sugar and 
phosphate (27). The force field for the TIS model includes 
base-stacking interactions which were taken from 
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empirical energy parameters (28). G»A pairs were given 
the same weight as A-U pairs. Replacing G»A parameters 
with those for G-C pairs resulted in negUgible differences 
to our simulation results (compare Figure 3 and 
Supplementary Figure S2). 

For thermodynamic analyses, multiple trajectories con- 
taining numerous folding and unfolding transitions were 
collected. All the simulations were performed using 
Langevin dynamics in the low friction limit (^ = 0.05 
^L~'). which enhances the rate of conformational 
samphng (29). The unit of time is = (ma^/Sh)'^^, in 
which the mass of the RNA beads is m = 0.1-0. 16 kg/niol, 
the average distance between beads is « = 4.6 A and the 
energy scale is £h ~1 kcal/mol. Profiles of the specific heat 
with respect to temperature, Cv(T), were obtained using 
Weighted Histogram Analysis Method (WHAM), and a 
peak in C/T) designates a melting transition temperature. 
The trajectories at the melting temperature were collected 
to generate the free energy profile with respect to the 
reaction coordinate, Q, the fraction of native contacts. 

RNA preparation 

Truncated P5abc RNA (6) containing 2-aminopurine 
(2AP) was chemically synthesized (Dharmacon) and 
purified by denaturing 10% polyacrylamide gel electro- 
phoresis. The RNA solution was desalted and 
concentrated by ultracentrifugation (Centricon Plus-20; 
Millipore) and stored at — 20°C in water. The typical 
yield after purification was 10-20%. 

UV thermal denaturation 

Thermal denaturation of P5abc was monitored at 
260 nm at a heating rate of 0.1°C/min on a Gary 400 spec- 
trophotometer. tP5abc RNA (2|iM) was dissolved in 
10 mM Na-cacodylate (pH 7.5), 40 mM NaCl plus the 
desired concentration of MgCl2. Global Melt Fit 
provided by D. Draper (http://ded.chm.jhu.edu: 16080/ 
—draper/) was used to obtain enthalpies and melting tem- 
peratures, assuming two transitions and temperature- 
dependent basehnes. 

RNase Tl partial digestion 

5'-pP]-labeled tP5abc was prepared with T4 polynucleo- 
tide kinase and gel-purified as described previously (30). 
The RNA was equihbrated lOmin in CE buffer 
(10 mM cacodylate, pH 7.5; 0.5 mM EDTA, pH 8.5) 
with or without 2.5 mM MgCl2 before treatment with 
0.1 U RNase Tl. Products were resolved on a 12% 
sequencing gel. 

Equilibrium Mg^* titrations 

For equilibrium folding, 2 ml 0.5 \\M RNA in CE and 2 ml 
0.5 |rM RNA in CE plus MgCl2 were denatured 5min at 
50°C and degassed. These solutions were mixed to 
produce Mg^^ titrations and incubated 5min at 30° C 
after each addition before the emission spectra (3-5 
scans) were recorded with excitation at 310 nm (Perkin 
Elmer LS50B). Excitation and emission slit widths were 
15 and 10 nm, respectively. 



The fluorescence intensity _F at 370 nm versus Mg con- 
centration, C, was fit to 

F = au+(«N+/^NC - (1) 

in which Cm is the midpoint of the folding transition, n 
is the Hill coefficient, is the fluorescence intensity at 
zero Mg"^ and a hnear baseline for the folded RNA is 
given by aN and j8n- The AG of folding was calculated 

from AG = RT^m^ ln(3+2V2), in which ^, = 
Cmax(4//<^C)max/^Catl/27w/.\' is a mcasurc of the 
cooperativity (31). Unfolding by dilution of Mg^"^ 
produced similar results (data not shown). 

Urea denaturation 

At equilibrium, 0.5 |iM RNA in CE plus IniM MgCl2 
was mixed with the same solution also containing 9M 
urea to produce the desired final urea concentration. 
The fluorescence data were fit to the linear binding 
model (32,33): 

(«N+/^Ni3)+(«u+/iu£))exp(-ff7(An - ^)/RT) 

l+exp(-77<Z)m-i')/RT) ' 

in which D is the concentration of urea, Z),n is the concen- 
tration of denaturant at the midpoint of the unfolding 
transition and m is a constant that reflects the change in 
folding free energy with denaturant concentration. Linear 
basehnes (a + pD) corrected for a shght change in the 
fluorescence of the folded state with urea, and a larger 
increase in 2AP fluorescence of the extended state with 
urea as stacking between 2AP and the adjacent 
bases was disrupted. The data were normalized as in 
Equation (1). 

Folding kinetics 

Folding kinetics was measured on an Apphed 
Photophysics SX.18MV stopped-flow spectrometer with 
a dead time of 1.3 ms. The excitation wavelength was 
310nm, and the emission signal was measured using a 
360 nm wide band cut-off filter. Folding reactions were 
performed by mixing tP5abc RNA prepared as mentioned 
earlier (1 |.iM in CE) with an equal volume of CE plus 
2 X MgCl2. Unfolding reactions were performed by 
pre-folding the RNA in CE plus 1 mM MgCl2 and 
diluting it 26-fold with buffer to reach the desired final 
MgCl2 concentration. Unfolding rates in urea were 
obtained by mixing 1 volume folded RNA (CE, 1 mM 
MgCl2) with 10 volumes urea in CE, 1 mM MgCl2. 
Folding rates were obtained by diluting RNA in 5.5 M 
urea 10:1 with buffer. Individual kinetic traces were fit 
to a single exponential rate equation and the resulting 
rate constants averaged, except for folding reactions of 
G174A, C166U RNA which were fit to a double exponen- 
tial equation. 

The observed rate Arobs as a function of Mg"^ concen- 
tration C was fit by a three-state model is ^ / ^ A^, in 
which the forward rate was given by /Cobs,f = kfC'\ where r 
is the apparent number of ions linked to the slow step of 
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Figure 2. Thermal stability of tPSabc RNA. (A-D) The first derivative of absorbance at 260 nm with temperature was fit to two transitions. All 
experiments were in 40 mM NaCl, lOmM Na-cacodylate (pH 7.5). Symbols: dotted fines, data; solid fines, first transition; long dashed lines, second 
transition; short dashed lines, baseline. (E and F) Simulated thermal denaturation profiles showing the specific heat C,, versus temperature for the 
folded (E) and extended (F) structures of tPSabc. 



the forward folding reaction (7). The unfolding pathway 
(A^ to /) was assumed to depend on the stability of N. 
Expressions for the forward and reverse paths were 
combined yielding the following equation: 




in which kf and /Cu are the folding and unfolding rate con- 
stants, K2 is the midpoint and m is the slope of the equi- 
librium between / and TV. 



P values 

The values of p for the Mg^^ titrations were calculated 
from j6 = n^/riii, in which «h = dln.Sr/dlnC was taken 
from the equihbrium folding curves and 
= dlnfcobs/dlnC was taken from the folding kinetics 
over the Mg^^ concentration ranges in which the folding 
and unfolding arms were hnear. 

RESULTS 

Thermal denaturation 

For these studies, we used the 56 nt tP5abc, which lacks 
five base pairs in the center of P5b and three nucleotides at 



the top of P5a (6) (Figure 1) yet forms the correct tertiary 
structure around the three-helix junction (12,17). UV 
thermal denaturation experiments in the absence of 
Mg'^"^ showed that the extended secondary structure 
unfolds in two transitions (Figure 2A) with midpoints at 
50 and 68°C (Supplementary Table SI). The melting 
profile in the absence of Mg"^^ was altered very little by 
the mutation A186U, which disrupts tertiary interactions 
between the A-rich bulge and P5c but does not change the 
predicted secondary structure (Figure 2B). In 2.5 mM 
Mg^"^, the secondary structure opened at a higher as 
expected (Figure 2C and D). Moreover, the 'wild type' 
tP5abc exhibited an additional low temperature unfolding 
transition not visible for the A186U mutant, which we 
assigned to denaturation of the tertiary structure. 

RNase Tl probing 

We next probed the tP5abc secondary structure by RNase 
Tl partial digestion at 25 and 50° C to verify that Mg^"^ 
induces the expected shift in P5c base pairing (Figure IC). 
At 25°C, which is below the first thermal transition, Gs in 
P5a, P5b and at the junction between P5b and P5c (G163 
and G164) were protected from cleavage, indicating that 
these helices are stable at 25°C in the extended structure 
(lane 3, Figure IC). As expected, these residues were 
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Figure 3. Free energy landscapes of tPSabc from simulation. Free energy profiles were projected to the fraction of native contacts (Q) for each region 
of tPSabc as defined in Figure 1. The free energy profile is calculated at the melting teinperature in units of k^T. (A-D) Energy profiles for the 
extended RNA. (E and F) Energy profiles for the folded state, at the high temperature transition (E) or low temperature transition (F). 
(G) Comparison of secondary structures (Qi) with native tertiary contacts, as reported by the overall root mean square deviation (dashed line). 



cleaved more strongly at 50°C, the Tm of the first thermal 
unfolding transition (lane 5, Figure IC). 

When the tP5abc was folded in 2.5 mM MgCU, 
G174-G176 in P5c and G180-G181 in P5a' also became 
protected at 25°C, showing that these base pairs near the 
three-helix junction become more stable in the tertiary 
structure (lane 4, Figure IC). Although bases in P5c are 
predicted to be paired in the extended structure, this helix 
may breathe or rearrange more than expected in low salt 
conditions. However, Mg^^ strongly enhanced cleavage of 
G169 in the loop of P5c, confirming the expected shift 
from a sheared G-A pair in the tetraloop of the 
extended form to an unpaired and exposed G in the 
pentaloop of folded tPSabc (lanes 4 and 6, Figure IC). 

Simulations of structural transitions during 
thermal unfolding 

To interpret the unfolding pathway implied by the thermal 
denaturation curves, we performed coarse-grained mo- 
lecular dynamics simulations of tPSabc using the TIS 
model in which each nucleotide is represented by three 
beads (27). We performed two sets of simulations, in 
which the TIS Hamiltonian was derived from either the 
extended tPSabc structure or the folded tPSabc structure 
(Supplementary Figure SI). In simulations of folded 
tPSabc, the specific heat (C,,) profile had a well-defined 
shoulder (Figure 2E), indicating two distinct melting tran- 
sition temperatures as observed in our experiments 
(Figure 2C). In contrast, the profile for the extended 
tPSabc had one peak (Figure 2F), corresponding to a 
single melting transition very similar to that of the 
A186U mutant in the presence of Mg^^ (Figure 2D). 
The overall specific heat profiles were in good agreement 



with the measured profiles, with similar transition tem- 
peratures, relative peak heights and overall shape. 

To quantify the structural transitions associated with 
the melting of tPSabc, we calculated a series of 2D, 
temperature-dependent free energy profiles, F(g„ g,), in 
which Qj and Qj are the fractions of native contacts of 
the individual helices in the extended tPSabc or folded 
tPSabc (Figure 3). Although the single peak in the C,, 
profile indicated cooperative folding of extended tPSabc 
(Figure 2F), the free energy profiles revealed intermediates 
corresponding to the formation of individual helices. The 
2D landscapes showed that the internal helices PSb and 
PSc form independently of each other with similar free 
energy barriers, in agreement with the RNase Tl data, 
while PSb is prerequisite for the closing helix, PSa 
(Figure 3A and B). PSa and PSc in turn stabilize PSa', 
completing the extended secondary structure (Figure 3C 
and D). Because the extended form of PSc is much more 
stable than the native PSc helix (Figure lA and B), our 
simulations predict that it forms early in the folding while 
native PSc does not. 

To compare the folding mechanism of the extended struc- 
ture to that of the folded structure, we constructed free 
energy landscapes for folded tPSabc at each of the two 
melting temperatures. These landscapes showed that the 
high temperature transition corresponds to sequential un- 
folding of PSa then PSb (Figure 3E), while the lower melting 
temperature corresponds to the cooperative unfolding of 
PSc and PSa' (Figure 3F). Thus, under native conditions, 
PSc and PSa' adjacent to the A-rich bulge fold coopera- 
tively, rather than sequentially as in the extended state. 

We next compared formation of tertiary structure, rep- 
resented by the overall root mean square deviation, with 
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Figure 4. Mg ^ dependence of folding at 30°C. (A) Cooperative change in the folding equihbrium with Mg ^. The fluorescence intensity at 370 nm 
versus Mg^^ concentration in lOmM CE was fit to a cooperative model [Equation (1)]; see Supplementary Table S2 for fit parameters. (B) Observed 
folding kinetics versus Mg~^ concentration. Errors are the standard deviation of three experiments. The line represents the best fit to the model in 
Equation (3), with = 36s^', k; = 4s^' (Supplementary Table S3). Filled and open symbols represent kinetics from folding and unfolding reac- 
tions, respectively (see 'Materials and Methods' section). (C) The midpoint of the Mg^^ titration was 10-fold higher in 60 mM Na^ than in CE buffer 
alone because the Na+ ions compete with Mg^^ for access to the RNA (Cm = 1.8 mM; Supplementary Table S2). (D) Observed folding rate versus 
Mg^^ concentration in 60 mM NaCl, as in (B), with A'„ s; 40 s^' and kf = 3.6 s^' (Supplementary Table S3). 



secondary structure formation as measured by g, 
(Figure 3G), and found that the tertiary structure forms 
in the lower temperature transition, concurrently with P5c 
and P5a'. Overall, the experimental data and the simula- 
tions show that native tP5abc unfolds through at least one 
equihbrium intermediate, corresponding to the loss of 
tertiary interactions and P5c. Our finding that the native 
P5c helix forms (or unfolds) in the same thennal transition 
as the tertiary structure shows that these steps are strongly 
coupled as previously suggested (6) and in agreement with 
the experimental RNase Tl results (Figure IC). 

Folding PSabc in Mg^^ 

We next asked whether similar intermediates were 
populated when P5abc refolds in solution. To study the 
Mg^"^-dependent folding pathway of tPSabc experimen- 
tally, we introduced the fluorescent base analog 2AP at 
position 185 (cyan; Figure lA). U185 stacks with 
adjacent adenosines in the extended conformer but flips 
out in the native structure. When U185 was replaced by 
2AP, the fluorescence intensity of the labeled RNA 
increased upon titration with Mg^^, consistent with 
unstacking of this base as the RNA folds 
(Supplementary Figure S3 A). 

Several control experiments indicated that the structure 
of tPSabc was not significantly perturbed by the 2AP sub- 
stitution. First, the 2AP substitution made almost no dif- 
ference to the thermal denaturation profiles of the RNA 
(Supplementary Figure S3B). Second, the folding time 
(?i/2 = 180 ms in 0.5 mM MgClz at 30°C; Supplementary 



Figure S3C) was compatible with previous estimates of the 
folding kinetics from time-resolved hydroxyl radical foot- 
printing (34) and NMR (12). Finally, a similar Mg^+-de- 
pendent folding transition was observed when pyrene was 
covalently attached to the 2' OH of U 142 (Suppleinentary 
Figure S4). 

The stability of tP5abc was obtained from the increase 
in fluorescence intensity with Mg^"^ concentration at equi- 
librium (Figure 4A). The midpoint of the folding transi- 
tion (Cm) was 0.18 mM Mg^+ at 30°C (in 20mM 
Na-cacodylate), reasonably similar to previous results 
(17). From the steepness of the folding transition with 
respect to Mg^^ concentration (31), we obtained 
AGuN = —1-3 ± 0.5kcal/mol (Supplementary Table S2). 

Mg^* dependence of the folding kinetics 

To probe the folding pathway of tP5abc in solution, the 
folding and unfolding kinetics was measured by 
stopped-flow fluorescence in Mg^"'" concentrations 
spanning the equilibrium transition (Figure 4B). The 
observed unfolding rate became faster as the Mg^^ was 
reduced below 0.2 niM, which was expected because the 
equilibrium no longer favors the folded state (Figure 4B). 
Surprisingly, however, the observed folding rate did not 
increase as the Mg^^ was raised above 0.2 mM, although 
the native tertiary structure was becoming more stable. 
These results showed that the folding rate is determined 
by a slow step that does not depend on the stability 
of the native state or the recruitment of additional Mg^"'" 
ions to the RNA. Higher salt did not eliminate the 
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Mg^"^-independent slow step in the folding pathway, as 
similar results were obtained when the tPSabc folding 
kinetics were measured in the presence of 60 mM NaCl 
(Figure 4C and D and Supplementary Table S3). 



B 



CD 

o 

o 
w 

CD 

o 

13 



CD 



CO 



173G A174 

t f A 
A A G G G U 

§ I • I I • 

G U U C C G^ 

\ I ° 

168C U166 



Mg2 



173G A174 
,t t A 

a'^aggg u 
c . I I 

GyU C CG 

* i ° 

168C U166 




WT 
U168C 
A173G 



50 
30 



10 
8 
6 



1 


1 1 1 J 1 1 1 1 1 III 


1 M 1 1 

WT 






U168C - 






A173G 










1 ' ' ' 1 ' 1 






0.1 


1 




[Mg2+] (mM) 





Figure 6. Rearrangement of P5c is not rate-limiting. (A) Base substi- 
tutions in P5c predicted to destabilize the extended conformation (left) 
without destabilizing the native structure (right). (B) Equilibrium Mg^^ 
titrations of wild-type tP5abc (black, filled circle), U168C (red, filled 
rhombus), A173G (blue, filled triangle) and G174A, C166U (green). 
Data were fit to Equation (1), which accounts for changes in 
cooperativity as well as midpoint (Supplementary Table S2). 
(C) Folding kinetics in Mg^^ at 30°C. Folding rates were measured 
as in Figure 4B and fit to Equation (3) (Supplementary Table S3). 
Open symbols, unfolding; filled symbols, folding conditions. 



Structural reorganization during folding 

To estimate how many interactions must reorganize 
during folding, we used an Arrhenius plot to measure 
the activation energy for folding in 0.5 mM Mg^"^, which 
was 21.3kcal/mol (Supplementary Figure S5). This acti- 
vation energy is comparable to AH for opening five base 
pairs based on empirical free energy parameters (35). 
Next, we measured the folding equihbrium and 
kinetics of tPSabc in urea. In 1 mM Mg^"^, urea denatur- 
ation of native tP5abc at equilibrium had a slope of 
m = 0.52kcal/mol»M (Figure 5A). Based on studies of 
RNA duplexes and tRNA, this sensitivity to urea is com- 
parable to the loss of four to live base pairs or the 
exposure of 2900 A^ of buried surface (33). Strikingly, 
urea strongly increased the folding rate under both 
native and denaturing conditions (Figure 5B). From 
these results, we concluded that tP5abc refolds through 
one or more intermediates and that formation of the inter- 
mediate requires the disruption of interactions in the 



extended state. However, we show below that this 
reorganization is not due to opening of P5c. 

Switching of the P5c helix 

Previous studies showed that mutations which favor the 
extended conformation of P5c inhibit base pair switching 
and tertiary folding, while mutations that favor the native 
secondary structure of P5c stabilize the native state (17). 
We used mutations that favor the native form of P5c to 
determine whether P5c switching slows down refolding 
of tP5abc (Figure 6A). These included U168C, A173G 
and the double mutation C166U, G174A, all of 
which disrupt base pairs in the extended P5c helix 
while maintaining the native base pairs. The mutant 
RNAs were shghtly more stable than the parental RNA 
(Figure 6B and Supplementary Table S2), confirming a 
very small shift in their folding equihbria toward the 
native structure. 
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Figure 7. Model for refolding of tPSabc. The folding pathway from the 
extended state ensemble {£} to the folded state N involves at least one 
metastable intermediate ensemble, {/}, which requires partial reorgan- 
ization of interactions in E and a second step that represents the Mg"^- 
dependent transition state (|). The switch in P5c occurs in the second 
step of the folding reaction, together with the formation of tertiary 
interactions between P5c and the A-rich bulge. Black, native conditions 
(ImM Mg^+); red, midpoint (O.lSmM Mg^+); orange, unfolded condi- 
tions (0.05 niM Mg^*). The landscape is based on the equihbrium 
stability of native tPSabc and simulations of experimental folding 
kinetics (Berkeley Madonna) assuming an intrinsic reaction time 
To = 1 us. The barrier between E and I was poorly constrained 
(dashed black line). 



Although we expected the mutant tP5abc RNAs to fold 
faster because the native P5c helix is more preferred, 
they refolded at a similar rate as the parental tP5abc (3 s~' 
versus. 4s~' in ImM Mg"^; Supplementary Table S3). 
About 15% of the C166U, G174A mutant refolded at a 
much slower rate (0.64 s~'; Supplementary Figure S6). 
Like the parental tPSabc, the folding rates did not increase 
with Mg ^ (Figure 6C). Thus, mutations that favor the 
native secondary structure of P5c failed to overcome the 
rate-limiting step of folding. From these results, we 
concluded that the overall folding rate is determined by 
reorganization of other interactions within the extended 
P5abc structure and not the P5c switch. In the subsequent 
sections, we discuss how the results of these experiments 
and our simulations can be rationalized by coupled forma- 
tion of tertiary interactions and helix switching. 



DISCUSSION 

How secondary and tertiary interactions are coupled in 
RNA remains an important question, as many RNAs 
such as riboswitches and even some ribozymes use helix 
switching to drive regulation or substrate recognition 
(36,37). When P5abc RNA folds, its secondary structure 
rearranges to a base pairing scheme that is less energetically 
favorable but that positions certain bases to participate in 
tertiary interactions (6). Mg^^ ions stabilize the tertiary 
interactions both directly and indirectly, compensating 
for the energetic cost of reorganizing the base pairs and 



maintaining a net favorable folding free energy (17,18). 
Thus, remodeling of the RNA secondary structure is also 
coupled to the specific coordination of metal ions. 

Tertiary folding drives helix switching 

Remarkably, both the simulated equilibrium folding land- 
scapes and the folding kinetics of P5c mutants suggest that 
P5c switches to its native form just before or concurrently 
with formation of tertiary interactions in the A-rich bulge 
(Figure 7). Not only is the native P5c helix predicted to form 
in the same thermal transition as the tertiary structure 
(Figure 3) but mutations that destabilize the extended 
form of P5c (and thus favor the native P5c) failed to 
speed folding (Figure 6), as would be expected if this 
switch happens before the tertiary structure forms. The 
folding kinetics also suggests that P5c switching occurs 
after the slow step of folding, which itself may involve 
local tertiary interactions around the three-helix junction 
(see below). Both aspects of our model are consistent with 
the proposal that the tertiary interactions directly promote 
the secondary structure rearrangement (18): the first step 
creates an intermediate conducive to switching, while the 
second step couples switching to tertiary interactions 
between P5c and the A-rich bulge. Ion concentrations, tem- 
perature and mutations, which modulate the folding land- 
scapes of P5abc, can be used to control the extent to which 
tertiary interactions determine P5c switching. 

In our model, the energetically costly P5c switch occurs 
very late in the folding landscape (Figure 7). As shown in 
our simulations, the native P5c helix is much less stable 
than the extended form and thus has httle probabihty of 
forming in the extended state. Thus, if the P5c switch 
occurred early, the activation barrier to folding would 
be very large, and P5abc might never fold. Instead, we 
suggest earher secondary and tertiary folding steps 
couple the change in secondary structure to the formation 
of tertiary interactions in the A-rich bulge, guiding the 
RNA through a lower energy path to the native state. 

Interestingly, the double mutation G174A, C166U, 
which stabilizes the native P5c helix relative to the 
extended form, results in biphasic folding kinetics. This 
could be due to the presence of an alternative folding 
path in which P5c switches before the rate-limiting con- 
formational change. Such a possibility is observed in the 
simulated folding landscapes (Figure 3B). 

Folding involves a metastable intermediate 

The tPSabc refolding kinetics in Mg^^ were simulated by a 
two-step folding pathway involving at least one intermedi- 
ate (/), consistent with high thermal activation energy of 
Mg^"^-induced refolding and the acceleration of folding by 
urea (Figure 7). SmaU differences in the observed rate con- 
stants of folding and unfolding reactions near the Mg^"'" 
midpoint (Figure 4B) suggest these pathways sample dif- 
ferent transition state ensembles consistent with this model. 

If P5c switching is not rate-limiting, folding must 
involve another structural change leading to the / state. 
This step could involve proper base pairing of P5a, which 
our simulations show depends on prior formation of 
the internal P5b and P5c hehces in the extended state. 
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The / state may also involve organization of the 
three-helix junction, including the tandem G-A pairs in 
P5b which are critical to the native tertiary structure. 
Substitution of G163 with 2,6-diaminopurine prevented 
tP5abc RNA from folding (Eda Koculi unpublished 
data), consistent with its participation in an extensive 
tertiary hydrogen bonding network. We considered 
whether the / state and the transition from I to N arises 
from reorganization of the tertiary interactions after 
(rather than before) P5c switching. However, this model 
does not explain why folding is insensitive to Mg'^"^ con- 
centration and is inconsistent with our simulations and 
thermal melting results. 

An extended transition state 

The relative position of the transition state ensemble along 
the reaction coordinate can be evaluated by comparing the 
sensitivity of the transition state free energy to Mg^"^ (n^ ) 
to that of the overall reaction free energy («h)' which is 
expressed as /S = «V«h (38,39). In unfolding conditions, 
/Su = 0.7-1, consistent with a transition state ensemble 
similar to the unfolded state, while under folding condi- 
tions, /ip = 0, which is also consistent with a transition 
state close to the unfolded state (Supplementary Table 
S3). These observations suggest that the transition states 
for both folding and unfolding are more open or extended 
than the native state and do not involve strong, 
site-specific coordination of metal ions. By contrast, the 
hairpin ribozyme folds through a compact transition state 
that is stabilized by Mg^^ similar to the native state (40). 
Thus P5abc, which must find a low energy path for helix 
switching, appears to adopt a different folding strategy 
than RNAs which fold via quasi-parallel docking of two 
hehcal domains. 

Effect on ribozyme folding 

Time-resolved footprinting experiments showed that the 
P5abc subdomain is the fastest region of the 
Tetraliymena ribozyme to fold (15). Therefore, the P5abc 
region hkely follows a similar refolding and hehx 
switching pathway in its native context as it does in isola- 
tion. However, interactions with other ribozyme sequences 
modulate formation of tertiary interactions in P5abc, and 
thus by extension, the rate of P5c switching. First, 
although the P5abc domain is marginally stable in isola- 
tion, extensive tertiary interactions between P5abc and the 
core of the Tetraliymena ribozyme stabilize the P5abc 
subdomain in its native context (41). Second, folding of 
P5abc is retarded in the wild-type ribozyme by premature 
docking against helices P4 and P6, suggesting that 
non-specific docking inhibits reorganization of P5c (34). 
Third, circular permutation of the ribozyme so that P5abc 
is near the 5'-end of the RNA causes P5abc to fold more 
rapidly than it does in the wild-type ribozyme (42). 

Although P5abc folds more rapidly when it is decoupled 
from other domains of the ribozyme, the wild-type con- 
nections improve folding of the ribozyme core (42,43). 
Not only does the tertiary structure of P5abc favor the 
active form of the ribozyme (13,44) but P5c switching is 
required for the long-range P14 kissing loop interaction 



between L5c and L2, which stabilizes the ribozyme core 
(45). As P14 impedes reorganization of the misfolded 
ribozyme core (46,47), the P5c switch might improve 
folding overall by delaying formation of P14 (12,48). 
Our work suggests that the switch kinetics is itself sensitive 
to alignment of helices in the ribozyme core, perhaps 
providing feedback on the quahty of the fold. 

CONCLUSION 

In conclusion, our simulations and experiments show that 
reorganization of the P5abc RNA secondary structure is 
tightly coupled to the formation of tertiary interactions, 
which are in turn stabilized by specifically bound Mg^"^ 
ions. Our results strongly suggest that these steps are kin- 
etically hnked and possibly concomitant, illustrating how 
coupled formation of secondary and tertiary structure in 
RNA can open up new, lower energy pathways to the 
native structure. Although simultaneous secondary and 
tertiary folding violates the hierarchy often assumed for 
RNA folding, it is commonly seen among proteins. Our 
results further imply that non-native interactions in the 
extended state aid folding rather than impede it. Because 
helix switching and P5abc tertiary interactions occur in the 
same step, this refolding pathway may improve coordin- 
ation between the stabilizing P14 interaction and proper 
docking of helices in the ribozyme core. Elucidation of the 
details of the folding mechanism and the extent of alter- 
native folding pathways will require a combination of 
ensemble and single molecule experiments and 
well-designed simulations. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Tables 1-3 and Supplementary Figures 
1-6. 
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